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Abstract 

The cross section for ultra-high energy neutrino-nucleon scattering is very sensitive 
to the parton distributions at very small values of Bjorken x (x < 10 -4 ). We nu¬ 
merically investigate the effects of modifying the behavior of the gluon distribution 
function at very small x in the DGLAP evolution equation. We then use the Color 
Glass Condensate formalism to calculate the neutrino-nucleon cross section at ultra- 
high energies and compare the result with those based on modification of DGLAP 
evolution equation. 


1 Introduction 


Neutrinos of ultra-high energies (E u > 10‘GeV) have been a puzzle for some time. One 
of the prime questions is where they come from, especially those above the GZK limit JT] . 
Possible sources include decays of super massive particles (dark matter?), acceleration in 
active galactic nuclei, and supernovae explosions (2|. Another question of interest is the 
cross section for the scattering of these ultra-high energy neutrinos with nucleons. Here, 
part of the interest stems from the fact that if the cross section increases sufficiently rapidly, 
then the unitarity limit may be reached .3j. Another interest is what one can learn about 
the very small x parton distributions, since the energy dependence of the inclusive cross 
section is very sensitive to them. 

The cross sections for scattering of neutrinos on nucleons at ultra-high energies are 
dominated by the gluons in the nucleon while the contribution of sea quarks is suppressed 
by a s since they come from gluon splitting via g —> qq. For x < ICC 2 the gluon distribution 
function of a nucleon is known to grow fast (3j with increasing Q 2 (virtuality of the gauge 
boson exchanged) and decreasing x as ( 1/x) 13 , with beta less than 1. This implies that the 
structure functions, e.g., F 2 , in deep inelastic scattering will also increase, which would in 
turn mean a fast increase of the neutrino-nucleon total cross section. This fast growth of 
the total cross section can not continue indefinitely since it would violate unitarity (the 
Froissart bound). The parton (gluon) phase space density (number of partons per unit 
area and rapidity) is expected to be very high at very small Bjorken x which would lead 
to an overlap in transverse space and recombination of gluons which in turn could lead to 
saturation (a slow down of the growth of the structure functions) and the unitarization of 
the cross section. 

At very small x , the nucleon is a very dense system of gluons and can be described via 
the Color Glass Condensate formalism |5j which resums large logs of energy as well as the 
large gluon density effects. It reduces to the BFKL formalism |B j in the limit that the gluon 
density in a nucleon is small. The Color Glass Condensate is an all twist formalism and as 
such extends the domain of applicability of pQCD to high gluon density environments. 

In this work, we consider different approaches to calculating the neutrino-nucleon to¬ 
tal cross section at ultra-high energies. First, we show the results from standard pQCD 
(DGLAP) |7j approach as well as the results from a unified DGLAP/BFKL approach, 
available in the literature |8|. We then consider the neutrino-nucleon cross section using 
the Color Glass Condensate formalism and gluon saturation based approaches. This in¬ 
volves modeling the quark-anti quark dipole cross section which is the basic ingredient in 
the structure functions. We compare the resulting neutrino-nucleon cross sections from 
different approaches and comment on the possibility of using future neutrino observatories 
to constrain the ultra-high energy neutrino-nucleon cross sections. 
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2 Neutrino-Nucleon Total Cross Section 


2.1 Leading Twist pQCD 


In perturbative QCD (pQCD), the cross section for the neutrino nucleon cross section can 
be written as 


vN 
Gtotal 



dQ 2 


d 2 a uN 
dxdQ 2 ’ 


( 1 ) 


where the differential cross section is given in terms of the quark and anti-quark distribution 
functions 


d 2 a uN 
dxdQ 2 


Gp / M\y, Z 
n \Q 2 + M'Q Z 


2 


q{x, Q 2 ) + (1 - Q 2 /xs) 2 q(x, Q 2 ) 


( 2 ) 


Here Gp is the Fermi constant and M wz refer to the W or Z boson masses while s is the 
neutrino-nucleon center of mass energy. The total cross section is finite (unlike the photon 
exchange process) and is dominated by scales Q ~ Mw,z■ In what follows, we will restrict 
ourselves to charged current exchanges, but the extension of work to the case of neutral 
current is trivial and we expect our results for the charged current exchange to hold equally 
well for the neutral current exchange. 

In the standard Leading Twist (LT) pQCD approach, one parameterizes the x depen¬ 
dence of quark and anti-quark distribution functions q(x, Q 2 ), q(x, Q 2 ) at some initial scale 
Q o, typically taken to be of order of a GeV or so. The distribution functions are then given 
by DGLAP evolution equations at any other x and 0 > Qq. The parameterizations are fit 
to the available data on DIS, for example, at HERA. There are various parameterizations 
of parton distribution functions satisfying the DGLAP evolution equations, for example 
CTEQ, MRST and GRV which differ in the choice of initial conditions and the degree of 
sophistication. 

If the neutrino-nucleon center of mass energy is much higher than the exchanged mo¬ 
mentum scale such that a s In s / M'Q ~ 1, it is more appropriate to use the BFKL formalism 
which resums these large logs rather than the DGLAP formalism. It is also possible to 
combine the two approaches in a phenomenological way such that both DGLAP and BFKL 
resummations are included |L In Fig. m we show the results of a DGLAP based cal¬ 
culation of the neutrino-nucleon total cross section, via charged current exchange due to 
Gandhi et al. J], denoted GQRS, as well as a calculation due to Kutak et ah, denoted KK 
(unified), which uses a unified DGLAP and BFKL approach (shown here without gluon 
saturation effects). The cross section grows with the center of mass energy which can be 
parameterized in terms of the incident neutrino energy (in the range shown in Fig. [TJ 
as cr ~ (E u /lGeV) 0A02 . It can be shown that this increase in the cross section is due to 
the growth of the parton distribution functions with decreasing Bjorken x 0. While at 
the lowest energy the two results are identical, which shows small x effects resummed by 
BFKL are negligible, at higher neutrino energy the two results can differ by a factor of two 
or larger. This signifies the fact that it is essential to include the contribution of small x 
partons properly at ultra-high energies. 
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Figure 1: The neutrino-nucleon cross section in Leading Twist pQCD via charged current 
exchange m 

It is important to realize that the HERA data on DIS covers a limited kinematic region 
and that ultra-high energy neutrino-nucleon cross sections are dominated by gluons at very 
small x and high Q 2 where there is no data. In the standard approach, one extrapolates the 
solution of the DGLAP evolution equations for parton distribution functions to smaller x. 
as needed. However, this requires making assumptions (or rather educated guesses) about 
the behavior of the distribution functions at small x. As we will show below, making rather 
plausible assumptions about the behavior of the parton distribution functions at small x , 
leads to large variations of the cross section at ultra-high energies. 

2.2 Gluon Saturation 

At very small Bjorken x , the gluon distribution function is expected to saturate, which 
would lead to a slow down of the growth of the neutrino-nucleon total cross section with 
energy. This is accomplished in the Color Glass Condensate (CGC) formalism which is an 
effective theory of QCD at high energy. The differential neutrino-nucleon cross section can 
be written in terms of the structure functions F\ and Fo (F :i does not contribute at small 
x )> 


with 

F 2 = 

Fi = 


d 2 


a 


G 2 


dxdQ 2 2vr (1 + Q 2 /M 2 V ) 2 


[(1 -y)F 2 (x,Q 2 ) + y 2 xF 1 (x,Q 2 )\ 


(3) 


J Q dz J dr 2 a d (x,r t ){Az 2 (l - z) 2 Q 2 K^(ar t ) + a 2 [z 2 + (1 - z) 2 )K 2 (ar t )} 
2x1^ Jo dZ I dr t ad (' Xj rt )a 2 [z 2 + (1 - z) 2 ]K 2 (ar t ) (4) 
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where a 2 = z( 1 — z)Q 2 and K 0 and K\ are modified Bessel functions, r t is the size of the 
dipole and z is the fraction of the photon energy carried by the quark. The total cross 
section is the integral of © over x, from x m i n = Q 2 /s to 1 and over Q , where we choose 
Qmin to be 10 GeV. The total cross section does not receive any appreciable contribution 
from scales below Qmin- The essential ingredient in saturation based approaches is the 
dipole cross section which is the imaginary part of the forward scattering amplitude (hence 
the name dipole cross section) of a quark anti-quark dipole on the nucleon. The dipole 
cross section ad(x,r t ) satisfies the JIMWLK evolution equation jTUj which is the all twist 
generalization of the BFKL evolution equation. In practice since the JIMWLK evolution 
equation is a highly non-linear equation, it is easier to parameterize the dipole cross section, 
in analogy to parameterizations of the standard parton distribution functions. The param- 
eterizations of the dipole cross section are then used to calculate the structure functions in 
© and checked against available data in DIS mmmm- The Color Glass Condensate 
formalism has also been successfully applied to particle production data in dA collisions at 
RHIC [ TO1 fTO . (for a review see P2P- The dipole cross section depends sensitively on the 
value of the saturation scale Q s and its energy dependence. While the overall magnitude 
of the saturation scale can not be determined from CGC itself, its energy (x) dependence 
is computed from CGC itself ca and is in good agreement with the value extracted from 
HERA phenomenology which has been parameterized m as 

Q 2 {x) = (!GeK 2 )(3 x 1(T 4 /x)' 28 . (5) 


The value of the saturation scale Q s compared to Mw determines whether one is in 
the saturation region (Q s > Mw ), in the so called geometric scaling |TDJ region (Q s < 
M w < Q 2 /A QC d) or in the DGLAP region ( Q 2 /A.q C d < M w )- It is ideal to have a unified 
formalism which can address all three regions; however, such a formalism does not exist 
currently. One can either use the DGLAP evolution equation and modify it to include 
gluon saturation effects as in 0 or use the CGC formalism and add the contributions of 
the DGLAP region by using the standard pQCD expressions. We choose the later approach 
since we are mainly interested in the ultra-high energy neutrino cross sections where the 
main contribution to the cross section comes from the very small x region. To do this, we 
introduce a cutoff x 0 below which we use the CGC expressions (J2J) while for x > x 0 we use 
© where the quark and anti-quark distributions are taken from CTEQ parameterization. 

One of the earlier parameterizations of the dipole cross section is due to Bartels et ah 
[©j which has been used to fit the HERA data. It is given by 

Od(x, rt) = 0 o[1 - exp(7T 2 r 2 a s (fx 2 )xg(x, /x 2 )/(3cr 0 ))] (6) 


with /i 2 = .26 /r 2 + 0.52 and the gluon distribution function xg satisfies the DGLAP 
evolution equation. The overal constant <7o is the nucleon size and taken to be <7o = 23 mb. 
This parameterization includes higher twist effects but does not have the BFKL anomalous 
dimension. Another parameterization is due to Kharzeev et al. na and has been used to 
fit the RHIC data on deuteron-nucleus collisions [is) nsi ■ The dipole cross section in this 
parameterization is given by 


<7d(r t ,y) = cr 0 




[rlQl(y)Y {vn) 



m 
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where the saturation scale is given by Q s (y) = Qo ex P[-Ml/ — Vo)/‘A with y = In 1 /a; and 
y 0 = 0.6, A = 0.3. The anomalous dimension 7 is 



where 



(9) 


This parameterization has the advantage that, unlike the one in (Rj|). it has the BFKL 
anomalous dimension built in which seems to be essential in describing the forward rapidity 
deuteron-gold data at RHIC. Using these two parameterizations of the dipole cross section, 
we calculate the neutrino-nucleon total cross section. We assume that quark (anti-quark) 
distributions are known well for x > Xq and use © to calculate the cross section for 
x > Xq. For x < Xq, we use the saturation approach and calculate the cross section using 
© with the structure functions given by using the two different parameterizations of 
the dipole cross section given in ©©, denoted BGBK and KKT dipoles respectively. To 
check the sensitivity of our results to the choice of Xo, we try two different values of xo, first 
Xq = lO -4 and then xq = 1CT 6 . In case of BGBK dipoles, since gluon distribution function 
xg{x, y 2 ) is not known well below x < 1CT 5 , we consider three wildly different scenarios; (i) 
a continually growing distributions for x < 10 -5 , (ii) a flat distribution for x < 10 -5 , and 
(iii) a distribution which falls by one order of magnitude for every decade of decreasing x 
for x < 1CU 5 . A measurement of the neutrino-nucleon cross sections at ultra high energies 
would thus go a long way toward understanding the very small x parton distributions. 

3 Results and Discussion 

In Fig. © we show our results for the neutrino-nucleon total cross section (via charged 
current exchange) for different neutrino energies for the case where xq = 1CU 4 and BGBK 
denotes the Bartels et al. model of the dipole cross section given by © and KKT denotes 
the Kharzeev et ah parameterization given in ©• The subscript / refers to the case where 
the gluon distribution function xg(x,y> 2 ) in @ , taken from CTEQ 6 , keeps growing with 
x below x = 1CU 5 while II refers to the case where the gluon distribution function below 
x = 10 -5 is flat and finally, case III corresponds to the case where the gluon distribution 
function below x = 1CU 5 falls like a power. 

For neutrino energies less than 10 s GeV, the cross section does not receive significant 
contributions from the region where x < 10 -5 . This shows in Fig. © as the three cases 
I, II, III (Bartels et ah dipole, denoted BGBK, with different gluon behavior at small x) 
being almost identical for E v < 10 8 GeV while the cross section calculated using the KKT 
parameterization of the dipole profile starts out below the other dipole models until about 
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Figure 2: The neutrino-nucleon cross section with x$ = 10 4 , details are given in the text. 


neutrino energy of 10 8 — 10 9 GeV after which it passes the BGBK //,/// dipoles, due to 
the constancy or drop off of the BGBK gluon distribution function below x = 10 -5 . 



Figure 3: The same as in Fig. © but with xq = 10 6 


To see the sensitivity of our results to the choice of cutoff x 0l we show the neutrino- 
nucleon cross section in Fig. © with the cutoff x 0 now taken to be 10 6 . Again, for x > x 0 , 
we use the quark and anti-quark distribution functions in © to calculate the cross section 
while for the region x < Xo we use the saturation approach. While the BGBK / does not 
change as it must not, the case where we have the gluon distribution function falling off 
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(BGBK III) is severely affected, by as much as a factor of 4 at the highest energy shown. 
On the other hand, the cross section using the KKT parametrization is rather robust, a 
change in xq from 10 -4 to ICG 6 changes the cross section by about 20% at E v = 10 10 and 
10% at E u = 10 13 . Depending on x 0 , the cross section given by Gandhi et al. |3] is about 
1.65 — 2.0 times bigger than the KKT cross section at E v = 10 12 . It is clear that the 
assumptions made on the behavior of the gluon distribution function at very small x will 
determine the outcome of the calculated cross sections at high energy. 



Figure 4: Ratio of KKT and KK (screened) cross sections . 

To compare our results to other saturation motivated studies, we show the ratio of 
our results for the neutrino-nucleon cross section using the KKT paraneterization, denoted 
KKT and the results of (screened) Kutak and Kwiecinski [SJ, denoted KK for the two 
choices of the parameter xq in Fig. Q. Since the numerical integrations involved are quite 
time consuming, we have have taken rather large increments in the integration rountines 
which leads to about 10% error on the KKT cross sections. This is the origin of the error 
bars shown in the figure. For neutrino energies more than 10 8 — 10 9 , the two approaches 
are in exccelent agreement for Xo = ICG 4 and within 20% for x = ICG 6 . The agreement 
is rather remarkable since the KK approach involves solving a phenomenologically unified 
DGLAP/BFK equation with a non-linear term motivated by the saturation physics while 
the results denoted KKT are based on a parameterization of the dipole profile which is 
motivated by the RHIC data on deuteron-gold collisions DU- This is most likely due to the 
similar growth of the saturaion scale in both cases since this growth is measures at HERA. 
It is also calculated very rcliabely in the Color Glass Condensate formalism |I2| and is 
in excellent agreement with the measured value at HERA. The fact that the two rather 
different approaches give quite similar results for neutrino-nucleon cross section at high 













energies is very reassuring and gives us confidence that if and when the ultra high energy 
neutrino-nucleon cross sections are measured, one can have quite stringent constraints on 
saturation based calculations of the neutrino-nucleon cross section. 
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